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I study the effects of disorder on the superconductivity of a three-band model with repulsive 
interband pairing. Such a model can support several possible superconducting order parameters, 
including a complex time-reversal symmetry breaking (TRSB) state. Impurity scattering suppresses 
the critical temperature of all these states, but the complex state survives, and remains a part of 
the phase diagram of the model even in the presence of moderate amount of disorder. This means 
that the TRSB states could be experimentally accessible in multiband materials like iron pnictides 
and chalcogenides. 


I. INTRODUCTION 


The discoveries of MgB^ 1 and the iron-based 
superconductors 2 have demonstrated that multiband su¬ 
perconductivity is a phenomenon with far reaching fun¬ 
damental and applied implications. It provides a new 
route to high-temperature superconductivity, and for en¬ 
gineering superconducting materials with novel proper¬ 
ties. Research in this area has been very active, and, as 
a result, rich physics has been uncovered. One of the 
most exciting ideas that have been discussed is the possi¬ 
bility of time-reversal symmetry breaking (TRSB) states 
in multiband systems 3 16 . They appear as complex ad¬ 
mixture of distinct and competing superconducting order 
parameters, like 5-wave and d-wave. The three-band sys¬ 
tem with repulsive interband pairing is the prototypical 
model which exhibits intrinsic frustration: the presence 
of several superconducting states with the same symme¬ 
try, which are directly competing with each other. In this 
model there are three possible two-gap order parameters 
with sign change between the gaps of each possible pair 
(s± states). They cannot be easily reconciled with each 
other, and the complex TRSB state appears as a natu¬ 
ral compromise for a three-gap order parameter. This 
model is thought to be relevant for the iron-based su¬ 
perconductors, most of which are believed to be in some 
form of the s± state (see recent review^®^). The TRSB 
state can potentially be tuned by doping in some of these 
materials 13 . 

Most of the research on TRSB states was done in the 
clean limit, ignoring the effects of the (unavoidable) im¬ 
purities and defects. However, it has been long recog¬ 
nized that even nonmagnetic impurities can be detri¬ 
mental to multiband superconductivity®^, especially 
in the case of an order para meter with sign change (for 
example, d -wave or s± state)PIIIS 28 . Effects of disorder 
on the TRSB states have not been studied so far, and 
it remains unclear if these states are reasonably robust, 
or, on the contrary, fragile and likely to be absent in 
any real material. In this paper I address this question, 
and consider the simplest possible three-band model with 
repulsive pairing interaction, which has well-established 
complex state in the clean limit. Introducing disorder, 


I study its effects on the critical temperature line and 
the possible phases close to it. I demonstrate that T c 
is strongly suppressed by disorder. The complex state, 
however, is still possible, almost until superconductivity 
is completely destroyed. This surprising, at first sight, 
robustness of the TRSB state is, in fact, easy to under¬ 
stand. As mentioned, this state is a complex admixture 
of three order parameter - the real two-gap s± states. 
Even though suppressed, these real states survive mod¬ 
erate amount of disorder. As long as they are present, 
there is always a region on the phase diagram where they 
compete, and in it the complex state emerges as a com¬ 
promise that minimizes the free energy. Thus, the TRSB 
state, arising in the tree-band model, could be observ¬ 
able in real materials (in which some degree of disorder 
is aways present). In contrast, for cases in which one 
of the competing phases is more susceptible to destruc¬ 
tion by disorder than the others (say d-wave vs. s- wave), 
it seems likely that the TRSB state region can be com¬ 
pletely wiped out. 

It is interesting to note that disorder-induced in¬ 
terband scattering can continuously change the order 
parameter of a two-band superconductor from s± to 
5 ++ 29 ^H state. This transition can involve an in¬ 
termediate complex combination of these two order 
parameter^, and also leads to TRSB state. It is ef¬ 
fect entirely due to the presence of interband impurity 
scattering. Thus, in some cases disorder is a necessary 
ingredient for the TRSB states. Furthermore, similar 
role to that of disorder can be played by a boundary 
scattering^, or a proximity-coupled gap!®®] 


II. THREE-BAND SYSTEM WITH 
INTERBAND PAIRING 


I consider a system described by the following Hamil¬ 
tonian: 
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Here i and j are band indicies, running from 1 to 3. 
The bands are parabolic, with partial and total densities 
of states (DOS) Ah, 7V 2 , N 3 and N t = Ni + N 2 + N 3 
respectively. For simplicity I assume they are identical, 
he., N 1 =N 2 = N 3 = Nq. 

The first term in Eq. [l] is the kinetic energy of the 
electrons, while the second term is due to impurities, 
which introduce both intraband and interband scatter¬ 
ing. The scattering rates are parametrized by 7 ^. The 
interband terms (i ^ j) can be written as 7 = NqT, with 
T = riimpTr^u 1171 ?) 2 , where 77 ^ mp and u irnv are the impuri¬ 
ties’ concentration and potential respectively (assuming 
Uimp is the same for each (i, j) combination). Since point 
defects, such as atomic substitutions or vacancies, can 
scatter carriers with large momentum change, the intra¬ 
band and interband scattering rates are expected to be 
comparable 37 . 

Superconductivity is driven by the pairing potentials 
TifJ, which can scatter Cooper pairs within the same band 
(i = j), or between different bands (i ^ j ). Note that 
the Cooper pairs themselves are always composed of elec¬ 
trons from the same band. In general, there are six in¬ 
dependent pairing interactions. For simplicity, I ignore 
the inraband terms u\f. Furthermore, I reduce the re¬ 
maining three interband terms u\^ to two, by assuming 
that U 12 = u± 3 . Finally, the interaction potentials are 
converted to dimensionless constants by using the DOS: 

A12 = N 0 u%, A 13 = N 0 ull A 23 = N 0 u s 2 c 3 - 
A12 = Ai 3 = A, A 23 = 77; A, 77 > 0 , 

The last inequality is true for systems with only repulsive 
interband pairing interactions. All the coupling constants 
can be combined into a 3 x 3 symmetric matrix 

/ 0 A A\ 

A $ A 0 77 . (2) 

\A V 0 / 


In the clean limit the gap parameters obey three cou¬ 
pled nonlinear self-consistency equations: 
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where uj = 7 rT( 2 n + 1) is the Matsubara frequency, and 
co>o is a high-energy cut-off (e.g., the Debye frequency). 
When disorder is present, it leads to additional inter¬ 
band scattering and coupling of the gaps, and the self- 
consistency equations become considerably more compli¬ 
cated (see, e.g., Ref. 138]) . However, close to the critical 
temperature these equations can be simplified by expand¬ 
ing them in powers of |A^|. In the presence of impurities 
this can be done systematically, starting from the Usadel 
equation d 25 * 39 ! 

The resulting gap equations can be though of as de¬ 
rived from an effective multigap Ginzburg-Landau (GL) 
free energy. I n the clean case this free energy looks rel¬ 
atively simpl^S 41 1 When interband impurity scattering 
is added, theory becomes more complicated; for the two- 
gap case it was derived earlier 32 . With three gaps the 
free energy can be written as 

gl = Ti + Tij + Fijk + Tem , (4) 

where we have kept only terms up to quartic order in A. 
The intraband terms 7y look like three copies of standard 
single-gap GL theory, 


^(A i )=a i |A ,| 2 + ||A i | 4 , (5) 

but with coefficients modified by disorder. In contrast, in 
Tij there are many new terms appearing due to scattering 
from impurities. It can be written as: 

Aj) = ajj|Aj||Aj| cos^jj + bij\Ai\ 2 \Aj\ 2 
+c u | Aj |' 3 1 Aj | cos 4>ij + Cij | Aj | 2 1 Aj | 2 cos ( 2 <y). ( 6 ) 


I study the phase diagram of the model, assuming that 7 
is fixed, and varying A from 0 (only two bands are cou¬ 
pled) to A 77 (a pair of bands is much weaker coupled). 
Of special interest is the region around the degeneracy 
point A = 77 . 

In the superconducting state there are three gap pa¬ 
rameters Ai, A 2 and A 3 , describing the superconduct¬ 
ing condensate in each band. They can be written as 
A i = \Ai\e l(pi ; the phases <pi are gauge-dependent, but 
their differences faj = cpi — (pj are physically meaningful, 
gauge-invariant quantities. Note that the superconduct¬ 
ing state breaks only the overall 17(1) symmetry. 

The system with A given in Eq. [ 2 ] has obvious frus¬ 
tration: each pair of gaps (i,j) prefers to have = 7 r, 

due to the repulsive Ay pairing. However, since there are 
three such pairs, there is no way to satisfy this preference 
for all gaps simultaneously. Thus, the different possible 
solutions compete, and this competition is the source of 
the TRSB state. 


For proc esses that contribute to these terms see Fig. [l] 
and Fig. [3p. In the clean limit ( 7 ^ —>> 0) the standard 
GL function for multiband superconductor is recovered: 
Tij dij\Ai\\Aj\ cos <pij. Notice that Tij is a simple ex¬ 
tension of the two-band theory 32 '. In contrast, the Tyy 



FIG. 1. Diagrams of processes contributing to the coefficients 
b \2 (on the left), and C 12 (on the right) of Jy?. The solid lines 
are electrons, and the dashed line represent the interband 
impurity scattering. 
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term contains entirely new, intrinsically three-gap con¬ 
tributions: 

Tijk ( 7 ^/c) dijk | A^ | | Aj | | A/^ | COS ((foij (/> 2 ^)( 7 ) 

This term is solely due to impurities, and it disappears 
in the clean limit (see Fig. [ 2 |. 




FIG. 2. Diagrams of processes contributing to the coefficients 
C22 (on the left), and /213 (on the right). The meaning of the 
lines is the same as in Fig. ^ 

Tem combines the electromagnetic field energy, and 
the derivative terms that couple A i to the electromag¬ 
netic vector-potential. For the rest of this paper uniform 
order parameter is assumed, and no field, so I will ignore 
Tem^ 

It is important to emphasize that the microscopically 
derived mean-field equations are the main object of our 
theory, and the free energy given in Eq. [ 4 ] is secondary. 
Although it may seem unnecessary to distinguish them, 
there are at least two problems in taking Tql as a start¬ 
ing point. First, the region of validity of the multiband 
GL theory is a matter of ongoing debatd^MJMII] j n addi¬ 
tion, and even more seriously, when the interband pair¬ 
ing dominates, this simple GL theory becomes unstable. 
This is due to the presence of a passive superconduct¬ 
ing channel, which has to be treated properly 15 ! The 
mean-field equations, however, remain valid. 


III. AT THE T c LINE 


At T c only the linear terms in the self-consistency equa¬ 
tions are important, and the problem can be reduced to 
a simple matrix one: 

MA = -A, ( 8 ) 

where I have defined the column vector A (with transpose 
A T = (A 1; A 2 , A 3 )). M is a 3 x 3 symmetric matrix: 


2A/y A/2 T 2\/y A/2 T 2X//y 

A/2+ 2\/ 7 (A + 77)^7 77/2 + (A + 77)/^ 

A/ 2 + 2\/ 7 77/2 + (A + rj)Ij (A + rf)Ij 


Here I have introduced 

UJQ 


-« LU (J 
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■27 


and / 7 = Ii — I2. The matrix is written in a form that 
separates the effects of the interband impurity scatter¬ 
ing. In the clean limit we have I2 —>• Ii = 70 ln(2a;o/7rT), 
(with 70 being the Euler constant), and / 7 0 , thus re¬ 

covering the equations obtained earlieJ 8 . Also, note that 
the intraband scattering rate does not appear at all - 
this is a consequence of the Anderson theorenP^. The 



FIG. 3 . T c for the two-gap and the three-gap order parameters 
as a function of the coupling constant A/77 f° r three disorder 
strengths 7/T c o: 0 (solid line); 0.2 (dashed line); 0.4 (dashed- 
dotted line). The critical temperature of the two-gap state 
does not depend on A. T c 0 is the value of T c for the two-gap 
state without disorder. 

eigenvectors of M represent the possible order parame¬ 
ters that can condense, and the corresponding eigenval¬ 
ues determine their respective transition temperatures. 
There are two possible order parameters, as in the clean 
case, with the third eigenvector having unphysical tran¬ 
sition temperature (as already mentioned, the presence 
of this passive superconducting channel has important 
consequence^). Quite surprisingly, the form of A at T c 
does not depend on 7 at all , and we recover the clean 
limit order parameters exactly: 


A«oc 

" 0 " 

-1 

, A (2) oc 

? 7+\/ 77 2 +8A 2 
2A 

1 


1 


1 


The first eigenvector represents a two-gap solution, in 
which there is a sign change between the two bands that 
become superconducting; the third band stays normal. 
The other eigenvector has gaps on all three bands, with 
two of them having the same sign. These solutions repre¬ 
sent different compromises for the underlying frustrated 
system. 
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The critical temperatures do depend, of course, on 7 , 
and are given by: 

- = h , 1 = 0 + 2 A - vV+8A)I 7 + T? ~ ^ 2+8A2 / 2 , 
77 z 

So how do the impurities affect the critical temper¬ 
ature? We plot the T c lines for both superconducting 
states as function of A for several values of 7 (shown on 
Fig. §. As we can see, the two states’ critical tem¬ 
peratures cross at A = 77 for any disorder strength. For 
A smaller than 77 the two-gap solution is the leading in¬ 
stability (since the system prefers to keep the weaker- 
coupled band normal). For A > 77 the three-gap state 
condenses first (and the weaker-coupled pair of bands 
has relative phase 023 = 0)- So far, this is completely 
analogous to the clean case. The only effect of impurity 
scattering is to suppress T c of both states, and eventually 
a critical amount of disorder destroys the superconduc¬ 
tivity in the system. This is as expected: because we 
are considering only repulsive interactions, superconduc¬ 
tivity depends crucially on the sign change in the order 
parameter. Since the interband scattering tends to av¬ 
erage the gaps on the different bands, large enough 7 
averages the order parameter down to zero. Plot of T c 
vs. 7 is shown on Fig. [4j While for the two-gap state the 
T c suppression follows the universal Abrikosov-Gor’kov 
(AG) curvd^, the three-gap case appears more compli¬ 
cated, with clear deviation from the AG curve at very 
low critical temperatures and large A/77 ratios. However, 
when A is close to 77 the disorder suppression of the two 
critical temperatures is (almost) equal. 

At the point A = 77 the two-gap and the three gap 
solutions are degenerate, and any linear combination 
of them is a possible order parameter. In particu¬ 
lar, a complex state can be constructed^ (A^ 3 )) T = 
A 0 (l, e 2?rz / 3 , e -27 ™/ 3 ). It was demonstrated in the clean 
limit that below T c this state plays an important role in 
the region of maximal frustration of the system A ~ 77. 
As we will see in the next section, this remains true in 
the presence of disorder. 


IV. TRSB STATE 

It was shown in the preceding section that there are 
two possible order parameters that can condense. Which 
one is leading depends on the ratio A/77. Below T c the 
situation becomes considerably more complicated, and 
when A ~ 77 the system may prefer a complex admix¬ 
ture of the two real states. To study this possibility 
we write a general complex order parameter (A^ 3 ^) T = 
Aq(— ft, e l(p , (it incorporates the 2 H 3 symmetry 

of the problem). Note that ip = 023/2. By appropriate 
choice of gauge Ai is made real (which is always possi¬ 
ble). Unlike the other two order parameters, however, 
A( 3 ) is intrinsically complex (for k / 0). This signifies 
the spontaneous breaking of the time-reversal symmetry 
of the system by the order parameter. It is caused by 



FIG. 4 . T c for the three-gap order parameter as a function 
of the disorder strength for five values of A/77 (k 1 order of 
increasing T c o): 0.90 (blue); 0.95 (blue-purple); 1 (purple); 
1.05 (purple-red); 1.10 (red). The T c lines for the two-gap 
and three-gap states coincide for A = 77 (purple line). T c0 for 
the two-gap case is the same as on Fig. [ 3 ] 


chiral pair currents in momentum space, when a pair is 
scattered from one band to another in a specific pattern 
(for example 1 -» 2 3). This symmetry breaking re¬ 

quires a second phase transition, below the customary 
superconducting one (associated with the 17(1) symme- 
try). 

It is important to understand the difference between 
the s ++ + is± states in the two-bancP^ and the three- 
band systems^. In the former the complex state is driven 
by interband impurity scattering, which leads to a quartic 
cos(2term in the GL expansion (a term that vanishes 
in the clean limit In contrast, the three-band case 
has an intrinsic instability towards a TRSB state, due 
to the competition between the three quadratic cos (4>ij) 
terms (see, e.g., Ref. my This explains the different 
behavior of the two systems close to T c (compare the 
phase diagrams constructed in Refs. j 8 j and my 

The intrinsic instability of the three-band model can 
be used to simplify the calculation. Close to T c the most 
important terms in the self-consistency equations are lin¬ 
ear (i.e., the terms with coefficients cq and a^). Thus, 
we can concentrate on the effects of impurities on those 
terms, and neglect all impurity induced terms, which are 
small anyway. (This is equivalent to only keeping 
from and completely ignoring 

With this, the first self-consistency equation is: 


cos ip = 


T]K 

2A 


/1 + 2A/ 7 \ 
VI + 277 TJ 


( 10 ) 


The other two equations are more illuminating if written 
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in terms of the GL coefficients cq and a^-, where we have 
= 7?/(2A 2 ) — /2, Cl2 = <23 = 1 /( 2 A) — I2 1 CL12 = ^13 = 
— (1/A + 2J 7 ), and <223 = — (I/77 + 2 / 7 ). (Note that < 2 ^ 
are negative, which seems to directly contradict the idea 
of phase frustration. It is, however, just a sign that the 
naive GL theory is unstable, and more care is needed in 
order to construct a proper GL free energy in the case 
interband pairing dominates 1 ^.) The remaining mean- 
field equations are 

K (fli = ^A^K 3 , (11) 


and 

a 2 

—ai(/^ 2 + 2)+/^ 2 12 +cos(2(^)a23 = 26 i(/^ 4 +2)Aq. ( 12 ) 

2(223 

Alternatively, these equations can be obtained from 
dT GL jdp = 0 , dT GL jdK = 0 and dT GL jd A 0 = 0 . 

Analyzing these equations provides a lot of useful in¬ 
formation. First, note that they reduce correctly to the 
clean limit results^ for J 7 —>> 0 . From Eq. [To] we see 
that nonzero k leads to <p 7 ^ tt / 2 . Thus, if we start from 
two-gap state (k = 0 , ip = tt/2) the system can continu¬ 
ously transform only to the TRSB state, but not directly 
to the three-gap one. Also, note from Eq. [IT] that the 
two-gap solution with k = 0 is always possible. However, 
if the right side is positive the equation can have also 
non-trivial k, solution. Since CL 23 < 0 the second term 
in the brackets is negative. On the other hand, a\ is 
positive at high temperatures, and only changes sign and 
turns negative for finite temperature below T c . Thus, the 
expression in the brackets can turn negative at some T 
below the critical temperature, and the TRSB state with 
n 7 ^ 0 and <p 7 ^ 7 r /2 emerges. 



FIG. 5. |A 0 (T)| (blue), k(T) (red), and tp(T) (green) for 

A /77 — 0.99, with 7 /T c o = 0.3 (moderate disorder). Imme¬ 
diately below T c the order parameter is in the two-gap state 
(k = 0, ip = 7t/2), but transitions to TRSB at T « 0.9T C . 


To demonstrate this, I solved Eqs. pR| pTT| and 12 


nu¬ 


merically, on both sides of the degeneracy point A = 77 . 
Results are shown on Figs. [5] and [ 6 ] For A < 77 the 


two-gap state is the first to condense, and for A > 77 
the three-gap state is leading. The TRSB state appears 
below T c in both cases, despite the presence of disor¬ 
der in the system. Note that the disorder is relatively 
weak, and the T c suppression is small - only about 30 
percent in both cases. However, calculations show the 
that the TRSB state survives even in the case of much 
stronger disorder, when the critical temperature is only 
a small fraction of the clean limit T c , provided that both 
real order parameters still have non-zero T c . Eventually, 
disorder completely suppresses the weaker superconduct¬ 
ing channel, at which point TRSB state also disappears. 
However, as long as the two channels are nearly degen¬ 
erate (i.e., A « 77 ), this single-order-parameter region is 
small, and TRSB state survives as a part of the phase 
diagram. 


K, ip 



FIG. 6 . The same as on Fig[5j but for A = I.OI 77 . The real 
three-gap state condenses first, but there is second transition 
to TRSB state around 0.9T C . 


V. DISCUSSION AND CONCLUSION 

In this paper I considered the influence of disorder on 
a simple three-band model with repulsive interband pair¬ 
ing interactions. This model is known to support a com¬ 
plex TRSB superconducting state in the region where the 
competition between the real two-gap and three-gap or¬ 
der parameters is strongest. As expected, the presence of 
impurities tends to suppress the critical temperature of 
all states, and for some finite value of disorder supercon¬ 
ductivity is completely destroyed. However, the complex 
state remains possible in parts of the phase diagram even 
in the presence of moderate disorder, as long as there are 
two real superconducting channels. 

The considered model and the derived results have to 
be taken with some caution. The three-band system is 
quite simple, and does not represent well the situation 
in iron-based superconductors (which, at the very least, 
have more bands at the Fermi level). This simplicity is 
the reason behind some artificial results, like the pres¬ 
ence of ungapped band in the two-gap phase. This phase 
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would likely disappear in a more realistic, and less sym¬ 
metric cases, in agreement with the general argument in 
Ref. [5(1 Furthermore, the results are derived using ex¬ 
pansion of the mean-field equations, valid only close to 
the critical temperature. In addition, I have neglected 
some quartic terms, with the assumption that they will 
not change the structure of the phase diagram (at least 
close to T c ). In spite of the simplifications, these re¬ 
sults are clear indication that TRSB states, if present in 
the iron-based superconductors, can be relatively robust 
against disorder, and thus experimentally accessible. 

The discussion in this paper has been mostly concerned 
with the stability of the TRSB state. Once it is present, 
equally important and interesting questions are of the 
physical properties and the possible manifestations of 
this state. Its presence requires a superconductor-to- 
superconductor transition, which may be detected in bulk 
thermodynamic property like the specific heat. Once in¬ 
side the complex state, its time-reversal nature can lead 


to induced currents and spontaneous magnetization close 
to impurities and edges^ (in the bulk the breaking is a 
consequence of interband Josephson currents, which do 
not lead to observable magnetization). Collective modes 
of the TRSB state are also quite interesting, as they con¬ 
sist of mixed phase-amplitude fluctuations of different 
gaps (see, for example, Refs. H5H5DE3. Another inter¬ 
esting consequence of the discreet symmetry breaking is 
the possibility of domain walls, whose presence can be 
used as a signature of this stated 
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